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(57) ABSTRACT 

A phasemeter for estimating the phase of a signal. For multi- 
tone signals, multiple phase estimates may be provided. An 
embodiment includes components operating in the digital 
domain, where a sampled input signal is multiplied by cosine 
and sine terms to provide estimates of the inphase and quadra- 
ture components. The quadrature component provides an 
error signal that is provided to a feedback loop, the feedback 
loop providing a model phase that tends to track the phase of 
a tone in the input signal. The cosine and sine terms are 
generated from the model phase. The inphase and quadrature 
components are used to form a residual phase, which is added 
to the model phase to provide an estimate of the phase of the 
input signal. Other embodiments are described and claimed. 

9 Claims, 2 Drawing Sheets 
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PHASE MEASUREMENT DEVICE USING 
INPHASE AND QUADRATURE 
COMPONENTS FOR PHASE ESTIMATION 

BENEFIT OF PROVISIONAL APPLICATION 5 

This application claims the benefit of U.S. Provisional 
Application No. 60/832,787, filed Jul. 24, 2006, which is 
hereby incorporated by reference. 

to 

GOVERNMENT INTEREST 

The invention claimed herein was made in the performance 
of work under a NASA contract, and is subject to the provi- 
sions of Public Law 96-517 (35 USC 202) in which the 15 
Contractor has elected to retain title. 

FIELD 

The present invention relates to signal processing and elec- 20 
tronic systems, and more particularly, to phase detection of a 
signal. 

BACKGROUND 

25 

Estimating the phase of a signal is an important function 
that finds application in a number of areas, such as in com- 
munication and signal processing. One such application in 
which phase estimation accuracy is important is in experi- 
ments for detecting gravity waves, whereby laser signals 30 
received from one or more spacecrafts convey information 
regarding small changes in the relative distances between the 
spacecrafts. Such changes in the relative distances may indi- 
cate gravity waves. For such an application, a locally gener- 
ated reference laser signal “beats” with the received laser 35 
signal to provide a signal comprising a tone having a fre- 
quency equal to the difference of the frequencies of the 
received and locally generated laser signals. This beat tone 
may be referred to as a heterodyne interference frequency, 
and the beat signal may be referred to as a heterodyne inter- 40 
ference signal. Accurately measuring the phase of the hetero- 
dyne interference signal provides information regarding the 
change in distance between the spacecraft. 

BRIEF DESCRIPTION OF THE DRAWINGS 45 

FIG. 1 illustrates a phasemeter according to an embodi- 
ment of the present invention. 

FIG. 2 illustrates a phasemeter for multi-tone input signals 
according to an embodiment of the present invention. 50 

DESCRIPTION OF EMBODIMENTS 

In the description that follows, the scope of the term “some 
embodiments” is not to be so limited as to mean more than 55 
one embodiment, but rather, the scope may include one 
embodiment, more than one embodiment, or perhaps all 
embodiments. 

FIG. 1 illustrates a phasemeter according to an embodi- 
ment. An input analog signal is provided to input port 101. 60 
The input signal may comprise one or more tones. For the 
present discussion, assume that the input signal comprises a 
single tone. (Embodiments for multi-tone input signals will 
be discussed later.) Under the assumption of a single-tone 
input signal , the input signal may be viewed as having a phase 65 
(|). Output port 120 provides an estimate of this phase, denoted 
as (j). 


2 

Although not necessary to the description of the embodi- 
ments, it is pedagogically useful to consider the physical 
significance of the phase term. In practice, the phase § may be 
viewed as the total argument of a sine function, and is a 
function of time. For example, a single-tone input signal, 
excluding noise terms, may be written as A(t)sin((|)). Here, 
is referred to as the phase of the input signal. Note that for the 
pure, single tone A sin(cot+(|)) where (|) is a constant, the phase 
is (j)=cot+(|), so that the time derivative of the phase is 00 , the 
frequency (in radians). 

Continuing with the description of the embodiment repre- 
sented by FIG. 1, the input signal is sampled by analog-to- 
digital converter 118 to provide a discrete-time, discrete- 
valued input time series, so that further operations are 
performed in the digital domain. This input time series is 
multiplied by two multiplier time series to yield what may be 
viewed as estimates of the inphase (I) and quadrature (Q) 
components of the input signal, where in FIG. 1 the signal 
flow path providing the estimate of the inphase component is 
labeled with an “I” and the signal flow path providing the 
estimate of the quadrature component is labeled with a “Q”. 
These two multiplier time series are cos((|)^) and sin((|)^), 
where (^is a phase term (it is a time series) provided by phase 
generator 1 02. In practice, these two multiplier time series are 
not exactly cos((() Af ) and sin^^) because once the input ana- 
log signal has been digitized, operations may only be per- 
formed with finite precision. With this in mind, the multiplier 
terms cos^^) and sin^^) will still be used, but it should be 
appreciated that embodiments will only utilize approxima- 
tions to these multiplier terms. 

Functional unit 103 and multiplier 104 indicate that the 
input time series is multiplied by cos((|)^), which is then low 
pass filtered and downsampled by functional unit 106. Some 
embodiments may not employ downsampling (decimation). 
For some embodiments, the low pass filter may be a Bartlett 
filter. Similarly, functional unit 105 and multiplier 107 indi- 
cate that the input time series is multiplied by sin^^), which 
is then low pass filtered and downsampled by functional unit 
108. Functional unit 108 has the same transfer function as 
functional unit 106. For some embodiments, a single func- 
tional unit may be used in a multiplexed fashion to implement 
the functional units of 106 and 108. For some embodiments, 
the cosine and sine terms may be evaluated by a table loop-up, 
with the model phase as an index (key) to a table entry, so 
that functional units 103 and 105 may be implemented by 
lookups into a cosine table and sine table, respectively. 

The output of functional unit 108 may also be viewed as an 
error term because it is provided to a feedback loop compris- 
ing controller 110 and phase generator 102. Although not 
necessary to the description of the embodiments, it is peda- 
gogically useful to consider the physical significance of this 
error term by ignoring for the moment that the input analog 
signal has been digitized and by ignoring that operations are 
performed in the digital domain. With this simplification in 
mind, if one considers a continuous -time input signal sin((|)) 
and continuous -time modeling phase then the error term, 

which is the low-pass filtered quadrature component of the 
term sin^sin^^), is given by sin^-^)^-^, assuming 
that is small. Consequently, the error term may be 

viewed as the phase difference 

Returning back to a description of the embodiment of FIG. 
1, functional unit 110 provides a correction term, denoted in 
FIG. 1 as A 00 , in response to the error term. For some embodi- 
ments, functional unit 110 may be implemented as a table 
look-up. For some embodiments, the correction term may be 
evaluated by multiplying the error term by a scalar. This 
correction term is added to an initial phase rate estimate oo 0 , 
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provided by functional unit 112 , to provide a phase rate term 
co. (This phase rate term is not to be confused with the use of 
co in the pure tone A sin (cot+(|)) example discussed previ- 
ously.) 

Phase generator 102 generates the phase term in 5 
response to the phase rate term co. The phase term may be 
viewed as a time integration of the phase rate term co, where 
because the operations are performed in the digital domain, 
the integration may be viewed as a multiplication by a dis- 
crete time index. That is, if an integer m keeps track of the 10 
number of clock cycles for phase generator 1 02 , then for some 
embodiments (|)^=oon. 

Because of the feedback loop comprising functional units 
106, 108, 110, the (model) phase term tends to track the 
phase of the input signal. However, utilizing the inphase 15 
component (the output of functional unit 106) in addition to 
the phase term ^provides additional accuracy in estimating 
the phase of the input signal. The inphase and quadrature 
components are provided to arctan 114. As its name implies, 
arctan 114 evaluates the arctangent of the quotient obtained 20 
by dividing the quadrature component by the inphase com- 
ponent, i.e., arctan(Q/I). This result may be termed a residual 
phase, denoted as A4> in FIG. 1 . This residual phase is added 
to and the result is low pass filtered and downsampled by 
functional unit 116 to provide the estimate <ji. Some embodi- 25 
ments may not perform low pass filtering, downsampling, or 
both, so that for some embodiments, <|> = <lw+A(|). 

Although not necessary to the description of the embodi- 
ments, it is pedagogically useful to consider the physical 
significance of the residual phase A(|) by again ignoring for the 
moment that the input analog signal has been digitized and by 
ignoring that operations are performed in the digital domain. 
With this simplification in mind, if one considers again a 
continuous-time input signal sin( 4 >) and continuous-time 
modeling phase then 


Q 

I 


sin (0-0 m) 
cos (0-0 M ) 


= tan {<p-<p M ), 


40 


so that A(|)=(|)-(|)^, which is recognized as a residual phase 
term. Of course, these expressions are not exactly satisfied by 
an embodiment, but are presented merely to understand the 45 
physical significance of the residual phase term and are not 
meant to describe the embodiment. 

Returning now to a description of the embodiments, the use 
of decimation (downsampling) filters allows for more than 
one clock domain. For example, for some embodiments, ana- 50 
log-to-digital converter 118, and functional units 102 , 103, 
104, 105, 106, 107, and 108 may operate at the highest clock 
rate; and the other functional units may operate at a lower 
clock rate. A potential source of error may be due to aliasing 
of noise during downsampling. This source of error may be 55 
mitigated by using digital anti-aliasing filters. For some 
embodiments, such aliasing filters are phase stable and linear 
in their passband. 

For some applications, the signal band may be very small 
compared to the heterodyne frequency, so that noise near the 60 
sampling frequency and its harmonics may alias into the 
signal band. As an example, for gravitational-wave observa- 
tion experiments, the signal band may be from 1 mHz to 1 Hz, 
so that for decimating from 10 kHz to 1 00 Hz, only noise near 
10 kHz and its harmonics is expected to alias into the single 65 
band. For some embodiments for such applications, by plac- 
ing nulls of the low pass and downsampling filter (e.g., func- 


tional units 106 and 108) at the sampling rate and its harmon- 
ics, aliasing of noise into the small signal band may be 
mitigated. 

For some embodiments, the low pass and downsampling 
filters 106 and 108 may be Bartlett triangular filters, which 
may be implemented by convolving two top -hat functions. 
Because the low pass part is followed by a decimation (down- 
sampling) stage, it is possible to save computations by calcu- 
lating only the points that will be retained, reducing the num- 
ber of calculations needed by the ratio of the output rate to the 
input rate. 

Under certain conditions, such as when the quantization 
error of analog-to-digital converter 118 is smaller than the 
other noise on the input signal, and when the signal frequency 
divided by the sampling frequency is not a rational number, it 
may be shown that the root power spectral density phase 
noise, as a function of sampling rate, f s , and number of bits 

N, is 


$=(ji2 N fiT s ) 

However, because the signal frequency may sweep through 
frequencies which are rational divisors of the clock, and the 
amplitude may not be exactly between two bit transitions, the 
above expression may not hold. To help mitigate phase noise, 
a small amount of noise in the least significant bit may be 
added for some embodiments. To generate such introduced 
noise, a pseudo-random number generator based on a linear 
shift feedback register, and a triangular distribution, may be 
used in some embodiments to dither the least significant bit. 

For some embodiments, the initial phase rate estimate pro- 
vided by functional unit 112 may be based on counting the 
zero crossings of the input signal. For some embodiments, the 
initial phase estimate may be entered manually, and for some 
embodiments, the carrier signal’s zero-crossing frequency 
may be added to a modulation offset frequency. Some 
embodiments may also include an automatic gain control to 
multiply the signal up to maximize the number of bits used in 
the calculation of the inphase and quadrature components. 

In some applications, the input signal may comprise more 
than one tone. To provide phase estimates for multi-tone input 
signals, the embodiment of FIG. 1 may be replicated, but 
where each replication has a different initial phase rate esti- 
mate. FIG. 2 illustrates such an embodiment, where there are 
K replications, and each box 201, 202, and 20 K is a replica- 
tion with a different initial phase estimate. For some embodi- 
ments, the approach for choosing the different initial phase 
estimates may be as discussed above with respect to the 
single-tone input signal. For some embodiments, a Fourier 
spectrum of the input signal may be evaluated, and the various 
initial phase estimates may be based upon recognizing the 
different peaks in the frequency domain. Other embodiments 
may employ other techniques for providing the various initial 
phase estimates, whether for the multi-tone case or the single- 
tone case. The above-described techniques are merely spe- 
cific examples. 

FIG. 2 also illustrates a space application of an embodi- 
ment for measurement of gravity waves, where a received 
laser and a locally generated reference laser are provided as 
inputs to photodetector 210 , followed by low pass filter 212 . 
The resulting signal provided to analog-to-digital converter 
214 is a heterodyne interference signal (at the heterodyne 
frequency). 

The various components and functional units described 
and illustrated in the various embodiments may be realized in 
hardware, software, or firmware. For example, some or all 
functional units may be implemented in one or more Field 
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Programmable Gate Arrays (FPGA), one or more Application 
Specific Integrated Circuits (ASIC), or some combination 
thereof. For some embodiments, one or more different clock 
domains may be utilized. 

Although the subject matter has been described in lan- 5 
guage specific to structural features and methodological acts, 
it is to be understood that the subject matter defined in the 
appended claims is not necessarily limited to the specific 
features or acts described above. Rather, the specific features 
and acts described above are disclosed as example forms of to 
implementing the claims. Accordingly, various modifications 
may be made to the described embodiments without depart- 
ing from the scope of the invention as claimed below. 

Throughout the description of the embodiments, various 
mathematical relationships are used to describe relationships 1 5 
among one or more quantities. For example, a mathematical 
relationship or mathematical transformation may express a 
relationship by which a quantity is derived from one or more 
other quantities by way of various mathematical operations, 
such as addition, subtraction, multiplication, division, etc. Or, 20 
a mathematical relationship may indicate that a quantity is 
larger, smaller, or equal to another quantity. These relation- 
ships and transformations are in practice not satisfied exactly, 
and should therefore be interpreted as “designed for” rela- 
tionships and transformations . One of ordinary skill in the art 25 
may design various working embodiments to satisfy various 
mathematical relationships or transformations, but these rela- 
tionships or transformations can only be met within the tol- 
erances of the technology available to the practitioner. 

Accordingly, in the following claims, it is to be understood 30 
that claimed mathematical relationships or transformations 
can in practice only be met within the tolerances or precision 
of the technology available to the practitioner, and that the 
scope of the claimed subject matter includes those embodi- 
ments that substantially satisfy the mathematical relation- 35 
ships or transformations so claimed. 

What is claimed is: 

1. An apparatus comprising electronic components to: 
multiply an input time series by a sine series to provide a 
first output; 40 

low pass filter the first output to provide a quadrature 
component; 

provide a model phase in response to the quadrature com- 
ponent, where the sine series is based on the model 
phase; 45 

multiply the input time series by a cosine series, the cosine 
series based on the model phase, to provide a second 
output; 

low pass filter the second output to provide an inphase 
component; 50 

provide a residual phase based on the arctangent of the 
quadrature component divided by the inphase compo- 
nent; and 


6 

add the residual phase to the model phase. 

2. The apparatus as set forth in claim 1, further comprising 
electronic components to: 

provide a correction phase rate in response to the quadra- 
ture component; 

add the correction phase rate to an initial phase rate esti- 
mate to provide a phase rate; and 

generate the model phase in response to the phase rate. 

3. The apparatus as set forth in claim 2, the apparatus 
further comprising electronic components to downsample 
when performing the low pass filtering of the first output and 
when performing the low pass filtering of the second output. 

4. The apparatus as set forth in claim 1, the apparatus 
further comprising electronic components to downsample 
when performing the low pass filtering of the first output and 
when performing the low pass filtering of the second output. 

5. The apparatus as set forth in claim 1, further comprising 
a sine table and a cosine table, wherein the sine and cosine 
series are obtained from a table with the model phase as a key 
to the table. 

6. The apparatus as set forth in claim 1, wherein the elec- 
tronic components are implemented in a field programmable 
gate array. 

7. A method comprising: 

multiplying an input time series by a sine series to provide 
a first output; 

low pass filtering the first output to provide a quadrature 
component; 

providing a correction phase rate in response to the quadra- 
ture component; 

adding the correction phase rate to an initial phase rate 
estimate to provide a phase rate; 

generating the model phase in response to the phase rate, 
where the sine series is based on the model phase; 

multiplying the input time series by a cosine series, the 
cosine series based on the model phase, to provide a 
second output; 

low pass filtering the second output to provide an inphase 
component; 

providing a residual phase based on the arctangent of the 
quadrature component divided by the inphase compo- 
nent; and 

adding the residual phase to the model phase. 

8. The method as set forth in claim 7, wherein low pass 
filtering the first output comprises downsampling, and 
wherein low pass filtering the second output comprises down- 
sampling. 

9. The method as set forth in claim 7, further comprising 
looking up a sine table with the model phase as a key to 
provide the sine series, and looking up a cosine table with the 
model phase as a key to provide the cosine series. 





